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Abstract. Galaxy interactions are known to trigger starbursts. Young mas- 
sive star clusters formed in interacting galaxies and mergers may become young 
globular clusters. The ages of these clusters can provide clues about the timing 
of interaction-triggered events, and thus provide an important way to recon- 
struct the star formation history of merging galaxies. Numerical simulations of 
galaxy mergers can implement different star formation rules. For instance, star 
formation dependent on gas density or triggered by shocks, predicts significantly 
different star formation histories. To test the validity of these models, multi- 
object spectroscopy was used to map the ages of young star clusters throughout 
the bodies and tails of a series of galaxy mergers at different stages (Arp 256, 
NGC 7469, NGC 4676, Arp 299, IC 883 and NGC 2623). We found that the 
cumulative distribution of ages becomes shallower as the stage of merger ad- 
vances. This result suggests a trend of cluster ages as a function of merger 
stage. In NGC 4676 we found that two clusters have ages of about 170 Myr, 
suggesting that they likely formed during its first passage. Their locations in the 
tidal tails are consistent with the spatial distribution of star formation predicted 
by shock- induced models. When comparing the ages and spatial distribution of 
clusters in NGC 7252 to our model, we found that some clusters are likely to form 
during the prompt starburst at first passage, as predicted by simulations with 
shock-induced star formation. These simulations show that the shock-induced 
mechanism is an important trigger of star formation and that using the ages of 
clusters formed in the starbursts can effectively determine the star formation 
history of merging galaxies. 



1 Numerical Simulations of Galaxy Mergers With Star Formation 



Galaxy collisions provide a natural laboratory for probing how star formation 
is affected by major rearrangements in the structure and kinematics of galac- 
tic disks. Many observational studies have been devoted to investigating these 
phenomena and helped to e stablish the link between galaxy interact ion and 
induced star formation (e.g. Kennicutt 19981 ; Sanders &: Mirabel 19961 ). How- 
ever, the triggers of star formation in interacting galaxies are still not fully 
understood. Studies have suggested two mech anisms to describe the star for - 
mation enhancement — density-dependent (e.g. Schmidtl 1959; K ennicutt]|l998l ) 
and shock-induced (e.g. lJog &; Solomonlll992l ; IScoville et al. 19861 ) star formation 
rules. Numerical models implementing these rules suggest that simple density- 
depende nt rules cannot offe r a complete description of star formation in merging 
galaxies ( Mihos et al. 19931). and that the two rules predict significantly different 
star formation histories ( Barnes! |2004| ). 
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Figure 1. Bolometric luminosities as functions of time. Top: shock-induced 
and Bottom: density-dependent simulations. Colors indicate galaxy pairs col- 
liding with different geometries: Black as DIRECT encounter, re d as POLAR, 
green as INCLINED and blue as RETROGRADE (|Barnesll2002l ) . Dashed line 
represents the first passage at T = 0. 



Fig. Q] shows bolometric luminosities as functions of time for merger simu- 
lations using t he two star formation rules and a sample of encounter geometries 
(|Barnea [2002) . Bolometric luminosity is a good tracer of the star formation 
rate (SFR) since it comes largely from young massive stars. Different encounter 
geometries yield different star formation histories, but the choice of star for- 
mation rules is clearly a more important factor. In shock-induced simulations 
(top), a global burst is triggered by large-scale shocks during the first passage 
at T = 0; later bursts of star formation, concentrated within the central regions, 
occur with the second passage and merger at ~ 500 — 600 Myr. In contrast, 
density-dependent models (bottom) generally predict a rather gradual increase 
within the central regions of the galaxies following the first passage; only the 
low-inclination passage (RETROGRADE) shows a starburst before the galaxies 
fall back together and merge. 



2 Young Star Clusters In A Series of Merging Galaxies 

Violent interactions often trigger starbursts which lead to the formation of young 
massive star clusters. These are likely to become youn g globular clust ers (YGCs) 
if they are still gravitationally bound after ~ 40 Myr (Schweizer 19991 ). The ages 
of these YGCs can be interpreted to yield the timing of interaction-triggered 
events, providing a powerful way to reconstruct the star formation history of 
merging galaxies. IChien et al.l ( 20071 ) obtained spectra of 12 young clusters in 
NGC 4676 using LRIS on Keck. These spectra yielded reliable approximations 
for cluster age and metallicity. Among the ages obtained, two are ~ 170 Myr, 
which suggests that they likely formed during the first passage of NGC 4676 
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(|Barnedl200? ). These two objects are located in the tidal tails of the pair, 
which is consistent wit h the spatial distribution of star formation predicted by 
shock- induced models ( Barnes! 120041 ). 



We have also obtained ages of clusters in a series of merging galaxies, rang- 
ing from early stages (Arp 256, NGC 7469) through merging (Arp 299) to fully 
merged (NGC 2623, IC 883) systems. For example, Fig. [2] shows spectra of 
6 young clusters in the merged system IC 883. Based on our age results, we 
compare the age distribution of clusters in each galaxy (including NGC 4676) 
according to their stage of merger (Fig. [3]). We found more than 70% of the 
observed clusters have ages less than 10 Myr in the first two mergers, Arp 256 
and NGC 7469, indicating strong on-going star formation in these galaxies. The 
ages are distributed more evenly out to about 260 Myr in the last two merger 
remnants, IC 883 and NGC 2623, which may suggest that some of these clus- 
ters formed during the first or second passages. This result suggests a trend of 
cluster ages as a function of merger stages: the cumulative distribution of ages 
becomes shallower as the stage of mergers advances. These age distributions 
provide a crucial way to discriminate between the alternate star formation his- 
tories predicted by the two rules described in Sec. 1. Detailed analysis of ages 
and metallicities of young star clu sters in these galaxies will soon be published 
( Chien et al. in preparation! 120091 ). 





hi 1 1 i i 


£* 

III 1 




II 1 1 1 1 1 



f 



-o .a or- 







......... 







3500 4000 4500 5000 5500 



3500 4000 4500 5000 5500 



Figure 2. ACS/WFC F435W image of IC 883 and spectra of 6 young clus- 
ters obtained with LRIS on Keck. Spectra are plotted as relative flux vs. 
observed wavelength (A). Markers in the spectra are the Balmer series. 
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Figure 3. Observed cluster age distribution of each galaxy. For a given 
panel clusters are plotted according to their age, with the youngest aligned at 
the top and the oldest at the bottom. Red points represent clusters with their 
spectrum dominated by the Balmer emission lines; blue are those dominated 
by the Balmer absorption lines, and green are those that have composite 
Balmer features. 



3 Combining Observations with Simulations 



Using interactive software ( Barnes Hibbard 2009I ) to match dynamical models 



to the observed morphology and kinematics of mergers, we built new models of 



NGC 7252 with the two star formation rules described above ([Chien Sz Barnes submitted 
2009). In our models, this proto-elliptical galaxy formed by the merger of two 
similar gas-rich disk galaxies which fell together ~ 620 Myr ago. Fig. Ujshows the 
spatial distribution of stellar particles of different ages using the two rules. Al- 
though on-going star formation occurs in the central regions in both simulations, 
the shock-induced simulation predicts that the products of past interaction- 
induced star formation are also dispersed around the remnant and along the 
tails. 

In addition to comparing our simulations with the observed kinematics and 
morpholo gies, we performed a detail ed analysis of cluster ages in NGC 7252 
(Fig. E]). iSchweizer fc Seitzerl dlOQSl . hereafter SS98) found 6 young clusters 



which lie between 3 — 15 kpc (~ 10 — 35") from the center of NGC 7252 and 
have ages of ~ 400 — 600 Myr, indicating that they formed early in the recent 
merger; the ages of these clusters are plotted in Fig. [SJ From the top panel 
we see that the shock-induced simulation (gray) produces prompt burst at first 
passage ~ 620 Myr ago while SFR rises more gradually in the density-dependent 
simulation (black) and peaks ~ 100 Myr later. Some of the cluster ages have 
a wide uncertainty; both of our simulations successfully reproduce the range of 
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Figure 4. Best match of the simulations of NGC 7252 
jCh ien fc Barnes submitted! 120091 ). Left: density-dependent and Right: 
shock-induced simulations. Old stellar particles are shown in grayscale. Red 
points are stellar particles with ages < 100 Myr, green with 400 — 500 Myr 
and blue with 500 — 600 Myr. The same number of points are displayed in 
both images. 



cluster ages, although cluster S105 and W6, with ages of 580 and 600 Myr, are 
more consistent with the prompt starburst at first passage in the shock-induced 
simulation. 

However the spatial distribution of the observed clusters strongly discrimi- 
nates between our models. To compare with the locations of the SS98 clusters 
(10" < r < 35"), the bottom panel of Fig. [5] shows the age distribution of 
star particles, located within this annulus, from our simulations. The density- 
dependent simulation (black) produces a gradually declining distribution of ages, 
with almost all interaction-induced star formation within the very central re- 
gions; for example the predicted SFR shows a broad peak around ~ 500 Myr 
ago but only a small portion of star particles within this annulus have such ages. 
In contrast, the histogram from the shock- induced simulation (gray) shows many 
more star particles in this annulus and a sharp peak around the first passage 
~ 620 Myr, suggesting that star formation occurs in more dispersed regions 
away from the centers and that a high portion of star particles within this an- 
nulus formed during the starburst in the first passage. This result may explain 
the distribution of ages observed in SS98, which indicates that shocks can be an 
important trigger of formation of these clusters. 

In summary, we show that besides the established role of density-dependent 
mechanism in enhancing the star formation rate in interacting galaxies, the 
shock-induced mechanism is another important trigger of star formation and 
that using the ages of clusters formed in the starbursts can effectively pin down 
the timing of interaction-triggered events and determine the star formation his- 
tory of merging galaxies. 
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Figure 5. Comparison of age distribution of stellar populations. Top Panel: 
Global star formation history (in simulation units) of NGC 7252 shown from 
900 Myr ago to present (T = 0). Black line repres ents density-dependent and 
gray shock-induced simulations. Cluster ages from lSchweizer fc Seitzerl ([19981 ) 
are plotted as dots with their uncertainties. Cluster SlOl has an age upper 
limit of 10 Myr and cluster S114 has an age of ~ 1 Gyr. Note that clusters 
S105 and S114 have possible ages of ~ 200 and ~ 40 Myr respectively. Bottom 
Panel: Histograms of number of stellar particles formed in the simulation, 
located within 10" to 35" from the center, measured at present time. 
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